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1.  Intrg<tucti9n 


The  fracture  toughness  of  certain  ceramic  composites  containing  second 
phase  particles  which  undergo  a  stress -Induced  phase  transformation  Is 
known  to  be  higher  than  that  of  the  brittle  ceramic  matrix,  [1,2].  A 
number  of  recent  studies,  beginning  with  the  work  of  HcMeeklng  and  Evans [3] 
and  Budlansky,  Hutchinson  and  Lambropoulos [ 4 ] ,  have  been  aimed  at  providing 
a  continuum  mechanical  model  which  can  predict  the  observed  stress 
Intensity  factor  reduction  at  a  crack- tip  In  such  materials;  see  also 
[5-9]. 

In  this  paper  we  examine  a  much  simpler,  but  related,  problem.  Ve 
consider  the  spherically  symmetric  deformation  of  a  hollow  sphere  which  has 
a  traction- free  Inner  wall  and  a  prescribed  radial  displacement  S  at  Its 
outer  wall.  Here,  one  expects  to  observe  a  reduction  In  the  stress 
concentration  factor  at  the  cavity  In  the  case  when  the  sphere  Is  composed 
of  a  transforming  ceramic  material  as  when  compared  to  the  case  when  it  is 
composed  of  the  ceramic  matrix.  We  obtain  a  closed  form  analytical  solu¬ 
tion  to  this  problem,  and  use  It  to  Illustrate  certain  features  of  boundar¬ 
y-value  problems  for  such  materials . 

The  analysis  here  Is  carried  out  within  the  small -strain  theory  of 
nonlinear  elasticity,  and  utilizes  the  particular  constitutive  law  proposed 
by  Budlansl^  et  al[4]  for  the  special  case  of  "supercritical 
transformations"  ;  this  constitutive  law  describes  a  class  of  elastic 
materials  which  Is  homogeneous  and  isotropic,  and  which  has  a  linear 
response  In  shear  and  a  trl-llnear  dllatatlonal  response.  Certain 


chcorecical  issues  pertaining  to  such  materials  were  examined  in  Part  I  of 
this  study,  [10]. 


We  show  that,  for  all  sufficiently  small  values  of  the  prescribed 
displacement  S ,  the  boundary -value  problem  has  a  unique  solution  and  that 
it  corresponds  to  a  configuration  of  the  body  in  which  the  strain  field 
varies  continuously  (a  "fully  untransformed  configuration");  this  is  also 
true  for  all  large  enough  values  of  S  (in  which  case  the  body  is  in  a 
"fully  transformed  configuration").  On  the  other  hand,  for  a  certain 
intermediate  range  of  values  of  S  the  problem  has  an  Infinity  of  solutions 
and  these  describe  configurations  which  involve  a  phase  boundary  ("par¬ 
tially  transformed  configurations”).  The  strain  field  is  continuous  on 
either  side  of  the  phase  boundary  but  suffers  a  Jump  discontinuity  across 
it;  the  displacement  field  is  continuous  everywhere. 


In  view  of  this  massive  failure  of  uniqueness,  we  are  led  to  conclude 
that  the  theory,  as  formulated,  is  deficient.  Presumably,  this  deficiency 
is  constitutive,  and  moreover,  is  Intimately  related  to  the  presence  of  a 
phase  boundary.  Accordingly,  in  Section  5  we  supplement  the  theory  with  an 
additional  constitutive  law  (a  "kinetic  law")  which  pertains  (only)  to 
particles  located  on  the  phase  boundary.  Since  quasi-static  motions  which 
involve  moving  phase  boundaries  are  generally  dissipative  (even  in 
nominally  elastic  materials,  see  [11]),  it  is  possible  to  define  a  notion 
of  a  driving  traction  on  a  phase  boundary;  the  kinetic  law  relates  the 
driving  traction  to  the  velocity  of  the  moving  phase  boundary. 


..4-- 

Uhlle  the  kinetic  law  governs  the  evolution  of  a  phase  boundary  once  it 
has  been  initiated,  a  separate  criterion  is  needed  in  order  to  signal  the 
first  appearance  of  the  phase  boundary.  The  "initiation  criterion"  used 
here  is  that  a  phase  boundary  will  emerge  when  the  driving  traction  on  it 
reaches  a  certain  critical  value.  In  the  present  theory,  this  is  equiva¬ 
lent  to  an  initiation  criterion  based  on  a  critical  value  of  the  dilata¬ 
tion. 

It  is  found  that,  as  the  given  displacement  6  increases  monotonically 
during  a  quasi-static  motion,  the  hoop  stress  at  the  cavity  first 
increases,  then  decreases  discontinuously  as  the  phase  boundary  emerges 
from  the  cavity  wall,  next  increases  slowly  (or,  for  certain  special 
kinetic  laws,  remains  constant)  as  the  phase  boundary  propagates  outwards, 
and  finally  commences  to  increase  at  the  original  rate  once  the  body  has 
been  fully  transformed.  In  general,  the  response  is  rate -dependent  and 
dissipative,  though  for  two  special  kinetic  laws  it  is  rate -independent  in 
one  of  these  special  cases  the  response  is  dissipation- free,  while  in  the 
other  it  is  "plasticity-like". 


Consider  a  hollow  sphere  of  incemal  radius  a  and  external  radius  b. 
Suppose  that  the  outer  surface  of  the  sphere  is  subjected  to  a  radial  dis¬ 
placement  S  while  its  inner  surface  remains  free  of  traction.  The  resulting 
deformation  of  the  sphere  is  assumed  to  be  purely  radial  with  u(r)  denoting 
the  radial  component  of  displacement;  u  is  required  to  be  continuous  on 
a^r^,  and  for  some  se(a,b)  it  is  to  be  twice  continuously  differentiable 
on  (a,s)+(s,b).  If  u'  is  discontinuous  at  r-s,  we  refer  to  the  circle  r-s 
as  an  equilibrium  shock  or  phase  boundary.  The  spherical  components  of 
strain  associated  with  this  deformation  are 

tee"  I 

V  for  r»-s,  (2.1) 

<r«  -  -  «rd  “  J 

and  the  corresponding  dilatation  d(r)  is 

A<r)  -  u'<r)  +  2  u<r)/r  for  r»*s.  (2.2) 


Suppose  that  the  sphere  is  composed  of  an  isotropic  elastic  material 
whose  stress-strain  relation  is 

2m€,+  (o(A)  -  2mA/3)1^  (2.3) 

fi  (>0)  is  the  shear  modulus  of  the  material  and  ff(A)  is  a  constitutive 
function.  Various  properties  of  this  class  of  materials  were  examined  in 
Part  I.  Here  we  simply  note  that  the  stress  response  of  this  material  in 
simple  shear  is  linear,  while  o(A)  denotes  its  stress  response  function  in 
pure  dilatation.  From  (2.1) -(2. 3),  the  components  of  stress  in  the  sphere 
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are 

^rr  -“ffCACr))  +  (4/i/3)  {u' (r)  •  u(r)/r}, 

<^i$  "  <^44  -  ^(A(r))  -  (2/i/3){u' (r)  -  u(r)/r),  •  for  Xf*s.  (2.4) 

°zB  "  <^64  “  "  0- 

Equilibrium  requires 
da^^/dx  +  2(oj.y  -  agg)/x  -  0  for  rf<s. 

Equation  (2.5),  in  view  of  (2.4),  (2.2),  leads  to 

^S(A(r))  -  0  for  rp^s, 

dr 

where  Z  is  defined  by 

S(A)  -  a(A)  +4/[iA/3  for  -«»<A<«b.  (2.8) 

(It  can  be  readily  shown  from  (2.8),  (2.3)  that  S  may  be  interpreted  as  the 
stress  response  function  of  the  material  in  uni-axial  deformation.)  Inte¬ 
grating  (2.7)  leads  to  2(A(r))  -  cj^  for  a<r<s  and  2(A(r))  -  C2  for  s<r<b 
where  cj^  and  C2  are  constants.  However,  as  shown  in  Part  1  of  this  study 
(see  (3.16)  of  [10]),  displacement  and  traction  continuity  across  r-s 
requires  S(A(r))  to  be  continuous  and  so,  (2.7)  in  fact  leads  to 

S(A(r))  -  c  for  r»<s.  (2.9) 

The  cavity  problem  consists  of  finding  a  displacement  field  u(r)  which 
obeys  equations  (2.2),  (2.9),  the  boundary  conditions 


(2.5) 

(2.6) 


(2.7) 
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u(b)  -  6, 

-  ff(A(a))  +  (4/*/3)  {u*  (a)  -  u(a)/a)  -  0, 


and  the  displacement  continuity  requirement  u(s-)-u(s+), 


(2.10) 

(2.11) 


In  this  section  we  will  solve  the  cavity  problem  for  the  particular 
class  of  materials  characterized  by  the  dilatational  response  function 


■  /3A  for  -An  £  A  ^  An, 

^(A)  -  ■  /9A  +  <7t(^'^m)/(^*%)  fof  %  A  <  Ajj, 

.  /3A  +  aj  for  A  >  A^; 


Aq,  An  and  are  material  constants  such  that 


(3.1) 


0>O ,  Ag|>An>0 ,  a'p<0 , 

•^An  <  /3Am  +  O-J  <  0 
(Affl-An) (^+4m/3)  <  -vj. 


(3.2) 


A  graph  of  the  function  o(A)  is  shown  in  Figure  1.  Requirement  (3.2)2 
implies  that  a(-An)  £  ^(Am)  <  0.  The  significance  of  (3.2)3 
discussed  shortly.  The  tri-linear  dilatational  response  function  (3.1)  was 
used  by  Budiansky  et  al[4]  for  characterizing  the  response  of  certain 
transforming  ceramics.  The  only  difference  between  (3.1),  (3.2)  and  the 
choice  made  in  [4]  is  that  we  take  o(A^)  to  be  negative;  this  feature  is 


needed  in  our  analysis  In  order  to  allov  for  the  occurrence  of  permanent 
deformations.  From  (2.8),  the  uniaxial  deformation  response  function  Z 
associated  with  (3.1)  Is 

od  for  -Afi  <  A  < 

2(A)-  *  aA  +  ax(A-Aii)/(Ajj-AH)  for  A^  <  A  <  A^,  (3.3a) 

^  oiA  +  <7^  for  A  >  Agj, 

where  a  -  d  (3.3b) 

The  third  condition  in  (3.2)  ensures  that  2(A^)  >  2(Aq)  so  that  2' (A)  is 
negative  on  (A^^.Ag).  As  shown  in  Part  I  (see  discussion  following  (3.16) 
in  [10]),  this  condition  is  necessary  and  sufficient  for  the  material  to  be 
able  to  sustain  equilibrium  deformations  with  discontinuous  strains.  (In 
the  terminology  of  Budiansky  et  al[4]  when  this  condition  holds,  the  mate¬ 
rial  can  undergo  a  "supercritical  phase  transformation".)  The  graph  of  2(A) 
is  shown  in  Figure  2;  the  number  2o  ■  (2f^+2„)/2.  The  figure  has  been  drawn 
in  the  case  2„^,  though  this  is  not  assumed  in  the  analysis. 

Finally,  we  introduce  some  additional  notation  which  will  simplify 
some  of  the  formulae  that  we  will  encounter  in  Sections  4  and  5.  This 
notation  pertains  to  certain  special  points  on  the  stress -strain  curves 
shown  in  Figures  1  and  2.  We  emphasize  that  the  constitutive  law  (2.3), 
(3.1)  involves  only  the  5  material  constants  n,  d.  A^,,  Af^  and  ar^\  the 
quantities  which  follow  can  all  be  expressed  solely  in  terms  of  these  basic 
quantities ; 
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<7K  - 

<^0!  “  (''ir‘^m^/2  -  aj(.l-fi/a)/2. 

+  aj, 

^1  “  ^  +  <^T;/a, 

^ol  “  (^+%)/2  +  oj/2a. 

Observe  that  the  points  (Ani.ffnl).  o"  first 

branch  of  the  stress-strain  curve  <7-o(A)  (hence  the  subscript  1),  while  the 
points  (Ag.ag,),  (^03,0^3)  and  (d){3,<7)(3)  lie  on  the  third  branch  (hence  the 
subscript  3).  Moreover  (d|(,<7)()  is  a  local  maximum  of  this  curve  (hence  the 
subscript  M)  while  is  a  local  minimum  (hence  the  subscript  m)  . 

Note  also  that  the  three  straight  lines  which  Join  (d^^.a^)  to  (d({3,C({3), 
<doi.<^ol)  (^I'^^ml)  (dn.^m).  each  have  the  same  slope 

-4/i/3.  The  points  (doi,<^ol)  and  (4,53,003)  correspond  to  so-called 
Maxwell -states;  they  have  the  property  that  S(doi)-r(Ao3)-ro-(2::m+%)/2 . 

We  now  return  to  the  cavity  problem  and  first  consider  the  case  in 
which  the  displacement  field  is  smooth.  Suppose  that  the  dilatation  d(r) 
is  such  that  <  d(r)  <  djf  for  a<r^,  so  that  all  particles  in  the  body 
are  associated  with  the  first  branch  of  the  stress -strain  curve  (i.e.  all 
particles  are  "untransformed").  By  (2.2),  (2.9),  (3.3)  it  then  follows 
that 

u'(r)  +  2u(r)/r  -  c/q  for  a<r^.  (3.5) 

Integrating  (3.5)  and  enforcing  the  boundary  conditions  (2.10),  (2.11) 


‘^ml  "  ‘^m  •  . 

‘^M3  ”  ''M  , 

“  ^^ol  - 

(3.4) 

%  -  <»dj4, 

dM3  “  ^M  ■ 

do3  “  (dm+dM)/2  -  o-r/la.  ^ 


lead*  to  the  following  expression  for  the  displacement  field: 


u(r)  -  (ib^/q)  {xr/a^  +  (l-ic)/r2)  for  a^r^,  (3.6) 

where  we  have  set 

a  -  W3a  (<1),  q  -  1  +  ((bVa^)-!)*  (>1) .  (3.7) 

On  using  (3.6)  and  (2.2),  the  requirement  lA(r)|  <  Afj  yields  Isj  < 
where 

6^  -  (qAMa3)/(3xb2).  (3.8) 

Next,  suppose  that  the  dilatation  is  such  chat  A(r)  >  A^^  for  a<r^  so 
that  all  particles  in  the  body  are  associated  with  the  third  branch  of  Che 
stress-strain  curve  (i.e.  all  particles  are  "transformed").  By  (2.2), 

(2.9)  and  (3.3)  it  then  follows  that 

u'(r)  +  2u(r)/r  -  (c-<rj)/a  for  a<r^b.  (3.9) 

The  displacement  field  may  now  be  found  by  integrating  (3.9)  and  enforcing 
the  boundary  conditions  (2.10),  (2.11).  This  yields 

u(r)  -  (Sb^/q)  {xr/a^  +  (l-a)/r2)  -  (<7'p/3aq){r  -  b^/r^)  for  a<r<b, 

(3.10) 

where  k  and  q  were  defined  in  (3.7).  On  using  (3.10)  and  (2.2),  the 
requirement  A(r)  ^  Aj,  yields  6  >  5^,  where 

■  (qAna^)/(3xb2)  +  .  (3.11) 

Finally,  we  consider  deformations  that  involve  a  phase  boundary  at  r-s . 


Suppose  that  A(r)>Aj,  for  a<r<3  and  that  |A(r)l<A(j  for  s<r<b,  so  that  all 
particles  vithln  the  phase  boundary  are  associated  with  the  third  branch  of 
the  stress -strain  curve  while  the  particles  outside  the  phase  boundary  are 
associated  with  the  first  branch  of  the  stress-strain  curve  (l.e.  the  body 
is  in  a  "partially  transformed  configuration"  with  the  particles  within  the 
phase  boundary  transformed  and  those  outside  untransformed).  From  (2.2), 
(2.9)  and  (3.3)  it  follows  chat 


u'(r)  +  2u(r)/r  - 


(c-aj)/a 
.  c/a 


for  a<r<s , 
for  s<r<b. 


(3.12) 


Integrating  (3.12)  and  enforcing  the  boundary  conditions  (2.10),  (2.11)  and 
the  displacement  continuity  condition  u(s-)-u(s+)  leads  to 


C  (Sb^/q)  {(cr/a^  +  (l-(c)/r2)  -  (tT^r/3aq) 

.  { (/eb^/a^)  (1-s^/b^)  (1-a^/r^)  +  (1-s^/r^))  for  a<r<s  , 


u(r)  - 


(3.13) 


(Sb^/q)  {icr/a^  +  (l-K)/r^)  -  (jco'jr/Saq) 

.(l-bVr^Xl-sVa^) 


for  s<r<b . 


The  restrictions  on  A(r)  that  were  assumed  in  deriving  (3.13)  can  now  be 
written  by  using  (3.13),  (2.2)  as 

%  -  aT(s3-a3)/(3ab2)  >  S  >  5^,  +  ax(b3-s3)/(3ab2) ,  (3.14) 

where  and  5^  are  given  by  (3.8),  (3.11). 


While  there  are  other  cases  to  be  considered  (for  example  A(r)>Aj5  for 
s<r<b  and  |A(r)|^f,{  for  a<r<s)  the  three  cases  considered  above  are  the 


2l.« 


P 

most  important  ones.  Arguments  similar  to  those  used  in  [12]  can  be  used 
to  show  that  the  remaining  cases  cannot  occur  in  any  quasi -static  motion 
H  which  commences  from  a  fully  untransformed  state,  provided  the  motion  con¬ 

forms  with  the  second  law  of  thermodynamics  under  isothermal  conditions. 


r 

I 


In  order  to 


E3  and  E3^  of  the  (£,s) -plane: 


suits  we  consider  the  following  sets 


El  -  {(5.s)|  l5|  <  Sfi,  s-a),  (3.16) 

E3  -  {(«,s)|  S  Z  S^.  s-b),  (3.17) 

E31  -  {(6,s)(  -  <7<j<(s^-a^)/(3ab2)  >  6  >  6^  +  ai«(b^-s3)/(3ab2) ,  a<s<b). 

(3.18) 

These  domains  are  sketched  in  Figure  3  where  Sj  is  defined  by 

S-j.  -  -  .  (3.19) 


While  the  figure  has  been  drawn  for  the  case  our  analysis  is  not 

restricted  to  this  case.  It  follows  from  the  analysis  surrounding 
(3.5)-(3.8)  that  if  (6,s)  is  a  point  in  Ei,  then  u(r)  as  given  by  (3.6)  is 
a  solution  to  the  cavity  problem.  Similarly  if  (6,s)  is  respectively  in  £3 
or  E31,  a  solution  to  the  cavity  problem  is  given  by  (3.10)  or  (3.13).  It 
is  now  clear  that  the  cavity  problem,  as  formulated,  suffers  from  a  massive 
failure  of  uniqueness.  Observe  from  Figure  3,  that  if  the  prescribed  value 
of  the  displacement  S  is  sufficiently  small  (-{[f  <  S  <  ,  the  problem 

has  a  unique  solution  (and  that  it  coincides  with  a  "fully  untransformed" 
configuration).  Similarly  if  S  is  sufficiently  large  (5  >  •  ^be 

problem  again  has  a  unique  solution  (this  time  corresponding  to  a  "fully 
transformed"  configuration) .  On  the  other  hand  when  the  given  value  of  6 


1 


lies  in  the  intezmediate  range  ^  <  £(f)-Sx>  problem  has  an 

infinite  nvunber  of  solutions  (since  the  value  of  s  is  essentially  arbi¬ 
trary)  . 


Hagp-atr^s  th?  gavUy 

Our  primary  interest  in  this  paper  is  to  examine  the  relation  between 
the  hoop  stress  at  the  cavity  wall  ■  agg(a)  and  the  applied 
displacement  S.  In  the  case  of  an  untransformed  configuration,  one 
finds  from  (2.4)2,  (3.1)  and  (3.6)  that  a^.  is  given  by 

^c  - 

Similarly  for  a  fully  transformed  configuration,  (2.4)2,  (3.1)  and  (3.10) 
give 

Oc  -  (3om/2«m)  «  +  3KOxbV2qa3,  (4.2) 

while  for  a  partially  transformed  configuration 

Oc  -  (3ffjt/2«M)  S  +  (3«<rx/2qa3)(*b3  +  (l-/e)s3),  (4.3) 

by  (2.4)2,  (3.1)  and  (3.13). 

The  regions  E3  and  E3]^  in  the  (6, s) -plane  are  carried  by  the 


respective  mappings  (4.1),  (4.2)  and  (4.3)  onto  the  following  domains 
F3  and  F3X  of  the  (S  ,<7g) -plane: 


r 


Fl  -  l(«.C7c)| 

-  6/i(l-»e)b25/a^q, 

1*1  ^  *M  1 ■ 

(4.4) 

F3  -  l(«.ac)| 

Of.  -  6fi(l-K)h^S/a^q  + 

(3»co3.b3/2qa^)  ,  f  > 

(4.5) 

F31  -  {(«.<rc)| 

Og  -  2k(Ty/2  <  2<r^5/2Sfn 

<  -  3/cofb^/2qa^ , 

(4.6) 

3^01/2  <  <^0  ^  3ffj£3/2  ) 

Figure  4  displays  these  regions;  and  F3  are  parallel  straight  lines, 
while  F3]^  is  a  parallelogram.  The  lines  s-constant  In  £3]^  are  mapped  onto  a 
family  of  parallel  lines  In  F33.  Observe  that  part  of  F3  coincides  with  one 
of  the  boundaries  of  F33,  but  that  this  is  not  so  of  F^^.  The  quantities 
^m*  ^T  given  by  (3.8),  (3.11)  and  (3.19),  while  the  numbers  a^, 

<7o3,  <^m3>  ^^ol  given  by  (3.4).  While  Figure  4  has  been  drawn  for 

the  case  0^3  >  0,  this  is  not  assumed  In  the  analysis. 
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5.  Klns.tisa 

In  order  to  complete  the  analysis,  we  must  account  for  the  kinetics 

-  -h  - 

of  the  transformation.  Let  A,  A,  Z,  Z  denote 

A-A(s+),  A-A(s-).  Z  -  Z  -  Z(A)  -  Z(A) .  (5.1) 


Then,  the  driving  traction  (the  driving  force  per  unit  area)  on  the  phase 
boundary  Is  (see  equation  (4.13)  of  Part  I) 


+ 

A 


f 


Z(A)dA 


A 


Z(A)(A-A) 


(5.2) 


On  using  the  constitutive  law  (3.3),  this  simplifies  to 


( 


f  -  (-a-^/a)  {  S  -  (%<-Z:„)/2) 


(5.3) 


in  the  case  of  a  partially  transformed  configuration  characterized  by 
(3.13).  Here  and  Sjj  are  as  defined  previously  in  (3.4)  and  have  the 
meanings  shown  in  Figure  2 .  Note  that  f  vanishes  when  Z  -  Zq  where 

Zo  -  (Z„  +  Z„)/2:  (5.4) 

Zq  is  called  the  Maxwell  stress  and  has  the  property  that  the  two  hatched 

+ 

areas  in  Figure  2  are  equal.  Since  Zq^<2^,  the  greatest  and  least 
values  of  the  driving  traction  (5.3)  are 

fn  -  -a^(5^-Zn,)/2a  (>0) , 

fm  -  aT<2M*2^)/2a  «0)  . 

respectively. 

Now  consider  a  quasi -static  motion  of  the  body  on  a  time  interval 
[tQiCi].  The  lack  of  uniqueness  observed  previously  suggests  that  the 
theory,  as  formulated,  suffers  from  a  constitutive  deficiency.  A  kinetic 
law  is  a  supplementary  constitutive  relation:  it  applies  to  particles 
located  on  the  phase  boundary  and  relates  the  driving  traction  f  (and  pos¬ 
sibly  other  local  quantities  as  well)  to  the  velocity  s  of  the 
phase  boundary.  An  example  of  such  a  kinetic  law  is 

s(t)  -  V(f(t))  for  to<t<ti,  (5.7) 

where  V  is  a  constitutive  function;  V  is  defined  and  suitably  smooth  on  the 


(5.5) 

(5.6) 


Interval  In  order  to  be  consistent  with  the  second  law  of 

thermodynamics  vuider  isothermal  conditions,  V  must  obey  (see  (5.2)  of 
Part  I) 

V(f)f  >  0  for  fa^f<fM.  (5.8) 

Returning  to  the  cavity  problem,  suppose  that  at  every  instant  t  during 
the  time  interval  [tQ.t]^]  the  body  takes  on  a  partially  transformed 
configuration;  the  displacement  field  in  the  sphere  is  then  given  by  (3.13) 
with  s  and  £  replaced  by  s(t)  and  £(t).  In  this  event,  we  find  from 
(3.13),  (2.2),  (3.3)  and  (5.1)  that 

21  -  3a£b^a/qa^  +  (ao-p/q)  (s^*^  *  I)  •  (5.9) 

Substituting  (5.9)  into  (5.3)  expresses  the  driving  traction  in  terms  of  S 
and  s.  Combining  the  resulting  equation  with  (5.7)  yields  the  following 
first  order  differential  equation  for  s(t): 

s(t)  -  V^(-a'j./a)  [  3a£(t)b2a/qa^  +  (aa<j./<l) (®^(t)/a^  -  1)  -  (5.10) 

for  tg^t^t]^.  Given  the  displacement  history  £(t)  for  to<t<t2^  and  the 
initial  position  of  the  phase  boundary  s(tQ),  (5.10)  can,  in  principle,  be 
solved  uniquely  for  s(t).  The  displacement  field  during  the 
quasi-static  motion  is  now  given  (uniquely)  by  substituting  this  s(t)  and 
S(t)  into  (3.13).  The  corresponding  history  of  the  hoop  stress  at  the 
cavity  is  likewise  given  by  (4.3).  We  now  consider  three  specific 
examples  of  kinetic  laws. 


Let  (p  be  the  function  which  is  inverse  to  the  kinetic  function  V; 
the  kinetic  law  (5.7)  can  then  be  written  in  the  alternate  form 

f(t)  -  <p(s(t))  for  to<t<ti.  (5.11) 

Consider  the  particular  kinetic  law  characterized  by 

^(s)  -  0  for  -«o  <  s  <  «o  (5.12) 

which  is  sketched  in  Figure  5;  according  to  this  kinetic  law  the  driving 
traction  f  on  the  phase  boundary  must  vanish  at  all  instants  during  a 
quasi 'Static  motion. 

First  consider  a  motion  which,  at  every  Instant  in  [tQ.tj^],  is  associ¬ 
ated  with  a  partially  transformed  configuration.  Eqtiations  (5.11),  (5.12), 
(5.3),  (5.4)  and  (5,9)  then  lead  to  the  following  relation  between  5(t)  and 
s(t) : 

3a5(t)b2#e/qa^  +  (saj/q) (s^(t)/a3  -  1)  -  Zg  for  to<t<ti.  (5.13) 

Eliminating  s  between  (5.13)  and  (4.3),  and  then  using  (3.4)  gives 

ac(t)  -  3ao3/2  for  tQ^t^ti.  (5.14) 

According  to  (5.14),  during  the  quasi-static  motion,  the  point  (6 (t) ,aj,(t) ) 
moves  along  the  horizontal  line  BC  in  Figure  5.  Recall  that  while  the 
figure  has  been  drawn  for  the  case  0^3  >  0,  this  need  not  be. 

Suppose  next  that  the  prescribed  displacement  5(t)  is 


increased  monotonlcally  (and  continuously)  from  zero,  and  chat  at  the 
initial  instant  the  body  is  in  a  £ully  uncransformed  configuration.  The 
resulting  history  of  the  hoop  stress  is  then  as  shown  in  Figure  5:  as  the 
point  (j(t) ,0(.(c})  moves  along  OA  (Che  sphere  remains  untransformed  and) 
increases .  When  ( (^ (t) ,<7g(c) )  reaches  point  A,  Che  particle  at  the 
inner  wall  r-a  is  at  a  "Maxwell  state"  in  the  sense  chat  Che  dilatation 
A  (a)  Ag2.  (so  chat  if  a  phase  boundary  was  initiated  at  r»a  at  this 
instant,  the  driving  traction  on  it  would  vanish).  If  we  assume  that  a 
phase  boundary  is  in  fact  initiated  at  this  instant  at  r-a,  (5(c),aQ(c)) 
now  moves  from  A  to  B.  The  hoop  stress  thus  decreases  discontinuously  from 
Che  value  3<rg]^/2  Co  3^03/2.  As  £(C)  continues  Co  Increase,  Che  motion  is 
now  governed  by  the  kinetic  law  and  so,  during  this  stage,  (5(t),ac(t)) 
moves  along  BC,  remains  constant,  and  the  phase  boundary  moves  outwards. 
Eventually,  (6(t) ,Oj.(t))  reaches  point  C  (at  which  time  the  phase  boundary 
has  arrived  at  the  outer  wall  r-b)  and  then  commences  to  move  along  C0«. 

The  hoop  stress  then  begins  to  increase  once  more. 

If  f(t)  is  decreased  monotonlcally  from  its  value  at  0^,  (5(t).ac(t)) 
follows  the  path  O^CBAO.  The  response  is  thus  reversible,  history* 
independent  and  dissipation- free. 


5.2  History- dependent,  rate -independent  response 

As  a  second  example  consider  the  following  choice  for  the  inverse 
kinetic  function  fp  in  (5.11), 


{f((  for  s  >  0, 

(5.15) 

for  s  <  0, 

which  Is  skecched  in  Figure  6;  here  f^  and  f^  are  the  maximum  and  minimum 
possible  values  of  the  driving  traction  as  given  by  (5.5),  (5.6).  Accord¬ 
ing  to  this  kinetic  relation,  In  order  for  the  phase  boundary  to  move  out¬ 
wards  the  driving  traction  f  must  take  on  Its  largest  possible  value  f^^, 
while  If  it  is  to  move  Inwards  f  must  have  Its  smallest  possible  value  f^; 
If  f  takes  on  any  value  between  f^  and  f^,  the  phase  boundary  must  remain 
stationary  (even  though  the  other  field  qtiantltles  might  be  varying) .  Thus 
(5.11),  (5.15),  together  with  (5.3),  (5.5),  (5.6),  (5.9)  and 
(3.4)  yield 


■  >  0 

if  - 

3<7j^3/2  and  6  >  0 , 

• 

s  * 

<  0 

if  Oc  - 

and  6  <  0, 

(5.16) 

.  -  0 

otheirwlse . 

As  the  following  two  examples  show,  the  response  of  the  body  to  various 
prescribed  displacement  histories  6(t)  may  now  be  determined  using  (5.16). 

Consider  first  a  monotonlcally  Increasing  displacement  history  6(t) 
with  6(tQ)-0.  Suppose  further  that  the  Initial  configuration  is  a  fully 
untransformed  one.  The  associated  variation  of  the  hoop  stress  Is  then  as 
shown  in  Figure  6:  As  (S  (t)  ,<7^(t))  moves  along  OVP,  Increases . 


When  (6(t) ,a{,(t) )  reaches  point  P,  the  dilatation  at  the  particle  at  the 
Inner  wall  r>a  Is  (and  thus.  If  a  phase  boundary  Is  Initiated  at  r-a  at 
this  Instant,  the  driving  traction  on  It  would  be  f^) .  If  we  assume  that  a 
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phase  boiuidary  is  in  fact  initiated  at  r-a.  at  this  instant,  (j(t),a^(t)) 
goes  from  P  to  Q  and  the  cavity  hoop  stress  decreases  discontinuous ly  from 
the  value  3ajj/2  to  3ajj3/2.  The  kinetic  law  (5.16)  governs  the  next  stage  of 
the  motion,  and  accordingly,  (6(c) ,ff(.(t) )  proceeds  along  QRS,  the  hoop 
stress  remains  constant,  and  the  phase  boundary  propagates  outwards.  Even¬ 
tually  (£(t) ,<7^(t))  reaches  the  point  S  (at  which  time  the  entire  sphere  is 
completedly  transformed)  and  then  commences  to  move  up  S0«;  thus  begins 
to  Increase  again.  If  6(t)  is  decreased  from  its  value  at  0*  the  path 
followed  on  the  (6 ,0^) -plane  is,  according  to  (5.16),  O^STUVO. 

Suppose  next  that  in  Che  preceding  example  the  displacement  6(t)  was 
only  increased  until  (6 (t) ,(t^(C) )  reached  point  R,  and  that  thereafter  it 
is  monoConically  decreased  for  a  short  interval  of  time.  According  to 
(5.16/,  s(c)  must  remain  constant  during  this  period  and  therefore 
(6(t) ,<7g(t))  moves  down  along  the  line  RW.  (RW  is  parallel  QU;  recall  the 
discussion  following  (4.6)).  If  6(t)  is  increased  again  from  its  value  at 
W,  (6(C)  ,ff(,(C) )  follows  the  path  VRS0«,  so  that  the  phase  boundary 
continues  to  remain  stationairy  for  a  while  (WR)  but  then  resumes  its  out¬ 
ward  motion. 

The  response  of  the  sphere  according  to  the  kinetic  relation  (5.15)  is 
thus  seen  Co  be  "plasticity-like".  Note  chat  quasi-static  motions  of  the 
sphere  are  dissipative  at  all  instants  during  which  the  phase  boundary  is 
in  motion  (since  then  fs>0)  but  non-dissipative  when  the  phase  boundary  is 
stationary.  This  particular  kinetic  law  is  equivalent  to  the  "flow  rule" 
used  by  Budiansky  et  al  [4] . 


As  a  final  example,  consider  the  kinetic  function  V(f)  shown  in  Figure 
7;  V  increases  monotonically  on  V(0)-0,  V(f)-Ko  as  f-»fM.  and 

V(f)-»-«o  as  f-»fni.  Suppose  that  5(t)-Xt  where  \>0  is  the  (constant)  loading 
rate.  During  the  resulting  motion,  the  point  (S(t) moves  along  the 
curve  OABMNO.^  shown  schematically  in  Figure  7.  The  sphere  initially 
remains  untransformed  (OZA) .  When  A(a)  -  (point  A)  we  assume  that  a 
phase  boundary  is  initiated  at  r-a;  the  driving  traction  on  this  phase 
boundary  at  the  instant  of  initiation  is  zero.  As  in  all  cases,  the  kinetic 
law  is  now  operative  and  governs  the  evolution  of  the  phase  boundary.  The 
equation  of  the  curve  BMN  is  found  by  solving  the  differential  eqxiation 
(5.10)  with  £(c)-At  subject  to  the  initial  condition  s»a.  It  is  clear 
that,  in  general,  different  loading  rates  X  will  give  rise  to  different 
curves  BMN.  Whan  the  rate  at  which  V(f)-»±«o  is  sufficiently  large,  one  can 
show  that  the  curve  BMN  does  not  intersect  the  upper  horizontal  boundary  of 
the  parallelogram  F3]^.  (If  V  does  not  increase  fast  enough,  the  path  will 
intersect  the  upper  boundai^;  this  means  that  the  sphere  cannot  be  deformed 
beyond  this  point  of  intersection,  at  that  same  rate  of  loading.)  Eventu¬ 
ally,  the  phase  boundary  reaches  the  outer  wall  (point  N) .  Note  that  dur¬ 
ing  this  motion,  the  cavity  hoop  stress  first  increases,  then  decreases 
discontinuous ly  as  the  phase  transformation  is  initiated,  then  increases 
slowly*  (as  the  phase  boundary  propagates),  and  finally  (once  the  entire 
body  has  been  transformed)  increases  at  the  same  rate  as  during  the  initial 
stage.  Unloading  follows  the  path  0.»CXYZ0. 

*  In  view  of  the  admissibility  condition  (5.8),  one  can  show  that  no  matter 
what  the  kinetic  law,  the  slope  of  the  curve  BMN  can  nowhere  exceed  the 
slope  of  the  straight  lines  OA  and  N0«. 
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Aa  a  second  example  of  a  loading  history,  suppose  that  the  Initial  con¬ 
figuration  of  the  body  is  that  associated  with  any  point,  say  M,  in  F3j^. 
Suppose  further,  that  the  displacement  6(t.)  is  held  constant  thereafter. 
Since  the  driving  traction  f  on  the  phase  boundary  does  not  vanish  in  gen¬ 
eral,  the  phase  boundary  will  move  according  to  the  kinetic  law.  The 
motion  of  the  body  may  be  determined  by  first  finding  s(t)  by  solving  the 
differential  equation  (5.10)  with  S(t)  -  constant,  and  then  substituting 
the  result  into  (3.13).  The  path  followed  in  the  (6, o^) -plane  is  the 
vertical  line  through  M.  The  phase  boundary  eventually  comes  to  rest  when 
the  driving  traction  becomes  zero  which  happens  when  (5(t) ,aj.(t) )  reaches 
the  line  BC. 
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Figure  3.  Parameter  sets  E^,  E^,  E3J  in  (6, s) -plane. 
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Figure  5.  Dissipation-free  response 
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on  the  phase  boundary  to  its  velocity,  is  then  imposed,  leading  to  a  unique  response 
in  all  quasi-static  motions. 

As  6  increases  monotonically  during  a  quasi-static  motion,  the  hoop  stress  at  the 
cavity  first  increases,  then  decreases  discontinuously  as  the  phase  boundary  emerges 
from  the  cavity  wall,  next  increases  slowly  (or,  for  certain  special  kinetic  laws, 
remains  constant)  as  the  phase  boundary  propagates  outwards,  and  finally  commences 
to  increase  at  the  original  rate  once  the  body  has  been  fully  transformed. 


